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Abstract The possibility of non-axisimmctric (kink) instabilities of a toroidal 
field seated in the tachocline is much discussed in the literature. In this work, the 
basic properties of kink and quasi-interchange instabilities, produced by mixed 
toroidal and poloidal configuration, will be briefly reviewed. In particular it will 
be shown that the unstable modes are strongly localized near the Equator and 
not near the Poles as often claimed in the literature. Based on the results of 
recent numerical simulations, it is argued that a non-zero helicity can already be 
produced at a non-linear level. A mean-field solar dynamo is then constructed 
with a positive a-effect in the overshoot layer localized near the Equator and a 
meridional circulation with a deep return flow. Finally, the possibility that the 
solar cycle is driven by a afl dynamo generated by the negative subsurface shear 
in the supcrgranulation layer will also be discussed. 
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1. Introduction 



An essential ingredient of any dynamo model is the a-effect, where a is the 
transport coefficient of the closure relation for the turbulent electromagnetic 
force. Owing to its pseudoscalar nature this term represents a likely possibility 
to produce a dynamo action with an axially symmetric field. 

From the theoretical point of view it is expected that both the kinetic helicity 
(v • V x v) and the current helicity (b • V x b) produ ced by the velocity v and 
magn etic field b fluctuations contribute to this term ( Gruzinov and Diamond! . 
1 9941 ) , although the relative importance of these terms is still a subject of de- 
bate ( Silant'ev . 200Cll ) . On physical grounds it is conceivab le that in the bulk of 
the convection zone "cyclonic" turbulence (jParkerl . 1 1 9 5 5h could be an efficient 
mechanism to p roduce a turbulent dyn amo whose a-effect is dominated by the 



kinetic helicity ( Steenbeck et ai , 19661 ) 



The discovery, due to the interpretation of hclioseismic data (|Schoilfl99lh . 
that the stable stratified region below th e convection zone is characterized by 
the presence of a strong horizontal shear (jKosovichevl . 119961 ) has suggested the 
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possibility of an a-effect located j ust beneath the convectio n zone (jParkerl . 11993 s ) 
in the region called the tachocline (jSpiegel and Zahnl . ll992f ). This fact has opened 
the door to the possibility that the source of the a-effect can have a magnetic 
origin, whi ch can thus be at tributed to various pos sible MHD instabilities fFerriz- 
Mas et al. Il994l; iDikpati and Gilmanl. l200ll ; IChatteriee et all l2004h . More 
recently, iBonanno and Urpinl (|201ll . l2012h " have further clarified a few aspects 
of kink and quasi-interchangc instabilities in stably stratified plasma. While 
kink instabilities are generated from a pure toroidal field, quasi-interchangc 
instabilities are produced by a mixed combination of the poloidal and toroidal 
field and their spectrum can be rather different from simple kink waves. It has 
also been shown that kink modes in spherical geometry are more effective near 
the Equator, a result that support the possibility of producing a non-zero a- 
effect at low latitudes. In fact, although the presence of a poloidal field breaks 
th e symmetry creates a preferred helicitv in the turbulence flow. Del Sordo et 
al. ( 2012h recently argued that kink instabilities alone can produce a preferred 
chirality in the turbulent plasma due to a symmetry-breaking effect at non-linear 
level. 

In this article, after reviewing the basic properties of kink and quasi-kink in- 
stabilities in a stably stratified plasma, the possibility of explicitly constructing a 
solar dynamo with an a-effect localized at the bottom of the convective zone and 
at low latitudes will be discussed and compared to mean-field mo dels where the 
a-effe ct is instead localized near the surface layer, as proposed by IBrandenbur j 
fl2005h . 



2. Quasi-Interchange Instabilities in a Stably Stratified Plasma 



The stability of a stably stratified column of plasma in the ideal MHD limit is 
the central problem of most of the co ntrolled- fusion literature. In this context, 
the energy principle ( Bernstein . 19581 ) has extensi vely be e n used in the past to 
study the stability of poloid al or toroidal fields ( Tavler . 1973al fbh and also of 
mixed combinations of both ( Tavlerl . Il980l ) . 

In cylindrical geometry, it can be proven that the plasma is stable for all 
azimuthal and vertical wave numbers ( m and fc), if it is stable f or m = in the 
k —> limit, and for m = 1 for all k (jGoedbloed et all 12004 ) . On the other 
hand, to show that a generic configuration with a combination of axial field 
and non-homogenous azimuthal field is stable against the m = 1 mode (for all 
k) is not an easy task in general and one has to resort either to a variational 
approach or to a numerical investigation of the full eigenvalue problem in the 
complex plane. In this respect, the "normal mode" approach can be more useful 
in astrophysics, as it is often important to know the growth rate of the instability 
and th e properties of the spectrum of the unstable modes (jBonanno and Urpinl . 
l2008al lbh. In particular, generic combinations of axial and azimuthal fields are 
subject to a class of resonant MHD waves that can never be stabilized for any 
value of the ratio of poloidal and toroidal fields. The instability of these waves 
has a mixed character, being both current- and pressure-driven. In this case 
the most rapidly growing unstable modes are resonant, i.e. the wave vector 
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k = (m/s)eg + k z e z is perpendicular to the magnetic field, B • k = where k z 
is the wavevector in the axial direction, to is the azimuthal wavenumber, and 
s is the cylindrical radius. The length scale of this instability depends on the 
ratio of poloidal and azimuthal field components and it can be very short, while 
the width of the resonance turns out to be extremely narrow. For this reason its 
excitation in simulations can be problematic. 

It is interesting to have a qualitative understanding of the MHD spectrum 
for a simple cylindrical plasma equilibrium configuration consisting of a mixed 
configuration of an azimut hal and a constant axial field B z . As was shown in 
Bonanno and Urpinl (2011), for a generic disturbance of the form e ( crt - lfc J 2 - mi v) 



an approximate expression for the dimensionless growth rate [r = a/u)A<t>] being 
u)A<j> the Alfvcn frequency in the azimuthal direction, is given by 

, 2m 2 (a - 1) 

(!) 



to 2 + (p 2 + m 2 )e 2 ' 

where p is the dimensionless radial wavenumber, e = B z /B v , a = din B v /dhis, 
and s is the cylindrical radius. Moreover the dimensionless vertical wavenumber 
[q = k z s] is close to the resonance condition [/ = qe + m « 0] which implies 
that the total Alfvcn frequency is zero [cja = B • k/-y/47rp = 0] at the resonance. 
The instability is never suppressed for any finite value of s and the growth rate 
is a rapidly increasing function of m in particular V 2 w (1 + e 2 ) -1 in the limit 
to ^> p 2 . If a < 1, it is possible to show that 

r 2 ~f 2 \±^, (2) 

1 — a 

that implies instability if a > — 1. The profile with a < — 1 is stable in this 
approximation. Note that modes with q satisfying the resonance condition uja = 
(or / = 0) are marginally stable because T = for them, but T 2 > in a 
neighborhood of the resonance. Therefore, the dependence of V on q should have 
a two-peak structure for any to. As in the case a > 1, the instability occurs for 
any value of e. If a = 1, then we have 



2to / 4to 2 



q 2 y (to 2 + q 2 ) 2 



(3) 



where fx is a positive number of the order unity. In this case, the dependence 
T 2 (q) also has a two-peak structure because T = at the resonance but T 2 > 
in its neighborhood. The instability is always present for any finite value of 
s. This explicit solution shows that, if a > — 1, the instability always occurs 
for disturbances with q and m close to the condition of magnetic resonance, 
uja = 0. The axial field cannot suppress the instability which occurs even if B z 
is significantly greater than B v . 



3. Kink Instabilities below the Tachocline 

Can the instabilities described in the previous section be operative in the over- 
shoot layer of the Sun, or even below? To address this question we must consider 
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the problem in spherical geometry, including the stabilizing effect of gravity and 
the destabilizing effect of thermal diffusivity. In fact the stability of the spherical 
magnetic configurations has been studied in much less detail and even the overall 
stability properties of stellar radiation zones are rather unclear. ? studied the 
stability of a random initial field in the stellar radiative zone by direct numerical 
simulations, and it was found that the stable magnetic configurations generally 
have the form of tori with comparable poloidal and tor oidal field strengths. Th e 
possible relaxation mechanism was further discussed bv lDuez and Mathisl (|2010h . 
The stability of azimuthal fields near the rotation axis has also been studied by 
Snruitl (Il999h . The author used a heuristic approach to estimate the growth rate 



and criteria of instability. Unfortunately, many of these estimates and criteria 
are misleading because they do not apply in the main fraction of the volume 
of a radiation zone where the stability pro perties can be qual itatively different. 
The heuristic approach was criticized by IZahn et al\ ( 2007|). The stab i lity o f 
the toroidal field in rotating s tars h as been considered by Kitchatinovl (2008) 
and iKitchatinov and Rudiger ( 2008 ) who argued that the magnetic instability 
is essentially three-dimensional and determined the threshold field strength at 
which the instability sets. Estimating this threshold in the solar radiation zone, 

the authors impose an upper limit on the magnetic field sa 600 G. 

The problem has recently been investigated bv iBonanno and Urpinl ( 2012 ) 
where it was shown that the most unstable modes have low r adial wavelengths 
at variance with the claim of IKitchatinov and Rudigerl ( 20081 ). Moreover if the 
thermal conductivity is considered, no threshold field is needed to trigger the 
instability. The most interesting result is illustrated by the angular dependence of 
the growth rate, as can be seen in Figure 1 where it appears that the instability is 
eff ective mostly n ear the Equator, and not along the axis, as originally supposed 
bv iTavler ( 1973af ). In general the location of the most unstable latitude depends 
also on the geometry of the basic state and cannot be determined only by local 
instability criteria. 

Can this instability produce an a-effect? It is difficult to consistently compute 
the strength of the a-effect within the linear analysis alone, altho ugh a first 
estimation in this direction has been provided bv lRiidiger et al. ( 2011 ). However, 
recent investigations in cylindrical symmetry pointed out that it is possible to 
generate a non-zero helicity out of a kink in stability due to non-linear effects 
(jDel Sordo et all 120121 ; IBonanno etall l2012h . 



4. Flux- Transport Dynamo and Tachocline a-Effect 

The previous discussion put forward that a likely location of the a-effect gen- 
erated by current-driven instabilities is precisely in the overshoot layer with a 
strong concentration of the turbulence at low latitudes where the instability is 
more effective. Although only 3D numerical simulations can in principle deter- 
mine the full structure of the a-tensor, as a first step in building a realistic model 
we can consider a latitudinal dependence of the type cos 9 sin 2 9 to model the 
suppression of the a-effect near the poles. The other essential ingredient is the 
inclusion of the meridional circulation. In fact, in the presence of a low eddy 
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Figure 1 . The dimcnsionlcss growth rate of the fundamental radial mode as a function of the 
polar angle, 9, for t = 5, m = 1, 2, and d = 0.1 in the case of a neutral stratification. Note that 
I is the latitudinal "quantum" number of the perturbation. Sec Bonanno and Urpin ( 2012) for 
further details. 



diffusivity [rj t ] (as is likely to be in the overshoot region), the magnetic Reynolds 
number becomes very large and the dynamics of the mean-field flow becomes an 
essential ingredient of the dynamo process. In this regime the advection produced 
by the meridional circulation dominates the diffusio n of the magnetic field which 
is th e n transported by the meridional circulatio n ( Dikpati and Charbonneatil , 
19991: iDikpati and Gilmanl. l200ll: iKiiker et "oil 



Chatteriee et all 



Pino,|2008|). 



20011 : 



Bonanno et ali. 



20021 : 



20041: iBonanno et ali 120061: Guerrero and de Gouveia Dal 



The magnetic induction equation reads 



dB 

~dt 



= Vx(U x B + aB) - V x (? ?T V x B) 



(4) 



where rj^ is the turbulent diffusivity. Axisymmetry implies that relative to spher- 
ical coordinates the magnetic field (B) and the mean flow field (U), respectively, 
read 



B = B(r, 9, i)e + V x [A(r, 9, t)e^} 



U = u(r, 9) + r sin 9Vl(r, 0)e 

being A(r, 8, t) the vector potential. The meridional circulation [u(r, 8)] and 
differential rotation [f2(r, 9)] are the poloidal and toroidal components of the 
global velocity flow field [U] . In particular the poloidal and toroidal components 
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of Equation (Q| respectively determine 



8A 1. _,. A . VT d 2 (rA) m d (ld(sA) 



(5a) 



dQd(Asine) ldfld(sA) 1 d / <9(rB; 

r 9r \ T 9?* 



— + sp(u V)f— ) 

9< \spJ dr 86 r 86 dr 



t]t 8 /ld(sB)\ 18/ 8(rA)\ 8 / a d(Asint 
~~^"dd\! 86 ) ~ rdrV dr ) ~ 86 \s^6 W 



(5b) 



where s = r sin 8. The a-effect and the turbulent diffusivity are parametrized by 
means of 



V 



1 2 

: — a>o cos 6 sin 6 
Vc + ^(Vt - Vc) 



1 



-erf 
erf 



Ol 



d 

r — r, 



1 - erf 



x - a 2 



(6) 



where ao is the amplitude of the a-effect, x = t/Rq is the fractional radius, oi, 
a,2 and d define the location and the thickness of the turbulent layer, rj t is the 
eddy diffusivity, r\ c the magnetic diffusivity beneath the convection zone and d v 
represents the width of this transition. In this investigation, the values a\ — 0.67, 
<i2 = 0.72, and d n = 0.025 have been used. 

The components of the meridional circulation can be represented with the 



help of a stream function ^ (r, I 



sin' 



6 cos 6 ip(r) so that 



9* 



r 2 psmd 86 ' 



rp sin 6 dr 



(7) 



with the consequence that the condition V-(pu) = is automatically fulfilled. 
A strategy to constrain several properties of the meridional circulation is to 
assume the differential rotation profile 0(r, 6) as a given ingredient, and deduce 
an approximation for the function ip from the angular-momentum conservation 
along the azimuthal direction. An approximate expression for ijj is thus 



bpr 



20, 



eq JO 



(u r ue)d6 



(8) 



where fi eq is the equatorial angular velocity. In part icular, for the standard 



isotropic mixing- length theory, Equation ([5]) becomes ( Durnevl . 120001 ) 



5 P r I 2\ 



20 



(9) 



cq 



In principle it would be possibile to explicitly compute ip and ip' using the 
relation ([9]) knowing the convective velocities of the underlying stellar model. In 
practice this would be problematic, because the convective fluxes and their radial 
derivatives computed from standard mixing-length theory are discontinuous at 
the base of the convective zone. In a more realistic situation the presence of an 



overshoot layer implies that 



smoothly so that ug is continuous at the 
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Figure 2. The a-effect (solid line), turbulent diffusivity (dot-dashed line) and (minus) the 
function ip(r) (dot-dot-dot dashed line), arc depicted in the left panel. The meridional circula- 
tion in units of the maximum surface value at a latitude of 45° is instead depicted in the right 
panel. 



inner boundary. Nevertheless one can use the representation @ to determine the 
stagnation point where ip' = 0, which turns out to be around x = 0.8 in terms 
of the fractional radius. An explicit form of the function ip which incorporates 
the following features reads 



ip = C 



1 — cxp 



(a; - x b f 



(x - 1) x 2 



(10) 



where C is a normalization factor, Xb = 0.65 defines the penetration of the 
flow, a = 0.08 measures how rapidly (uf) decays to zero in the overshoot layer 
and the location of the stagnation point. The density profile is taken to be 

— — xq ) in which m is an index representing the the stratification of 
the underlying solar model, its value in the region of interest is approximately 1.5, 
and xq = 0.85, so that with these value the strength of the meridional c irculation 
at low latitudes is of the same order as t h e surf ace flow, as discussed bv lBonanno 
(j2012l ) ; see also iPipin and Kosovichevl (|201ll ) for a similar investigation. The 
radial profile of the a-effect, turbulent diffusivity, stream function and meridional 
circulation used in the calculation are depicted in Figure [2] 

The hclioscismic profile for the differential rotation is taken so that 



n(r,9) =n c + s 



(O s (0)-O c ), 



(11) 



where O c /27r = 432.8 nHz is the uniform angular velocity of the radiative core, 
Q s (0) = fi C q + «2 cos 2 8 + a,4 cos 4 8 is the latitudinal differential rotation at the 
surface and 6(x) = (1+erf (x))/2. In particular f2 eq /27r=460.7 nHz is the angular 
velocity at the Equator, a2/27r=— 62.9 nHz, and 04/277=— 67.13 nHz. In this 
calculation the angular velocity is normalized in terms of equatorial differential 
rotation f2 oq , r c = 0.71 and d c = 0.025. As usual, the dynamo equations can 
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Figure 3. Butterfly diagram for a solution with meridional circulation included. The isocon- 
tours for the toroidal field (upper panel) at the base of the convection zone and the radial field 
on the top boundary are shown. Full and dashed lines show the positive and negative levels 
respectively. The solution is obtained with C a = 2, Cn = 4.5 X 10 5 and C u = 400 which implies 
a period of about 24 years a turbulent diffusivity r]t = 3.1 X 10 11 cm 2 s~ 1 and a (poleward) 
surface flow of about 18 m s — 1 . It is interesting to notice that at low latitude the phase relation 
is consistent with the observations, namely B r B^ < 0. This type of solution also shows a strong 
preference for dipolar modes. In particular the critical solution for the symmetric mode has 
C a = 3.70, significantly greater than the antisymmetric mode. An essential ingredient in order 
to get this property correct is the presence of a strong meridional circulation at the bottom of 
the convective zone. In particular in this case the bottom flow is about 8.5 m s _1 . 



be made dimensionless by introducing the dynamo numbers Co = R^Qeq/r/T, 
C a = RQCta/r/T, C u = RqUi/ttt, C u = RqU /rjr, where w is the frequency of the 
dynamo wave, and U = ug(r = Rq,9 = 45°). This linear dynamo problem is 
solved with a finite-difference scheme for the radial dependence and a polynomial 
expansion for the angular dependence by imposing potential field as a boundary 
condition for the surface, and perfect conductor for the inner boun d ary. Further 



details of t he numerical appro ach can be found in iBonanno et all (|2002l . [2006) 



and also in Jouve et al. (2008) 



A reference solution can be then obtained in the region of high C u , and its 
main properties are the following: C a — 2, Co = 4.5 x 10 5 , and C u = 400 which 
implies a period of about 24 years, a turbulent diffusivity i] t = 3.1 x 10 11 cm 2 s _1 
and a (poleward) surface flow of about 18 ms" 1 . The basic properties of this 
solution arc depicted in Figure 3 where the butterfly diagram for the toroidal 
field at the bottom of the convection zone, and of the radial field at the surface 
are shown. It is interesting to notice that at low latitude the phase relation is 
consistent with the observations, namely B r B^, < 0. The other property of this 
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solution is that it shows a strong preference for dipolar modes. In particular 
the critical solution for the symmetric mode has C a = 3.70, significantly greater 
than the antisymmetric mode. On the contrary, had we considered the case of an 
a-effcct uniformly distributed throughout the whole of the convection zone, the 
result would have been C Q = 2.1 for both symmetric and antisymmetric mode 
with a difference of about 1% most probably due to numerics. The conclusion 
is that the parity of the solution is more likely to be a dipole if the a-effect 
is located at the bottom of the convection zone. In our simulations the precise 
value of the ratio 7? c /?7t was not crucial as wc obtained basically the dynamo 
solution for rj c /r} t = 0.005 and r] c /r]t = 0.05 although in all of the simulations 
presented in this work the value Tj c /Vt = 0.01 has been chosen. For instance 
the period obtained for rj c /i]t = 0.005 was 24.9 years, and for T) c /r)t = 0.05 was 
24.5 years. This is not surprising as the penetration of the meridional flow is 
very weak as showed in the right panel of Figure © and therefore the precise 
value of rj c cannot significantly change the type of dynamo action. In addition, 
changing the width of the a effect and of the turbulent diffusivity also did not 
lead to significant changes in the solution, although we expect that the spatial 
extension of the turbulent layer should be of the order of the tachocline width 
d n w 0.05 solar radii. 

However, the serious, unsatisfactory aspect of advection dominated dynamo 
is the fact that the strength of the return flow is largely unknown, and the eddy 
diffusivity is about one order of magnitude greater than would be expected on 
the basis of the standard mixing length theory 



5. Solar Dynamo from Subsurface Shear Instabilities 

In recent years the possibility that the dynamo operates in the subphotospheric 
layers of the Sun mostly driven by the negative gradient o f the angular velocity 
near the surface has been proposed (jBrandenbura . |2005j). It is interesting to 
see if a mean-field dynamo model can describe this case with the help of a more 
refined differential-rotation profile including the negative shear in the subsurface 
layers. In order to discuss t his issue a slightly m odified version of the analytical 
approximation presented in lDikpati et al\ ([20021) has been used. The radial profile 
of the a-effect and turbulent diffusivity are depicted in Figure In particular 
the rotation rate is taken constant in the radiative interior [O c ] and the tachocline 
is located at the same radius as in the model described in the previous session. It 
is assumed that has constant width w O.O5i?0, and at the top of the tachocline 
its rotation rate is given by 



fi(r cz , 6) = fi C z + o-i cos 2 9 + ci4 cos 4 



(12) 



where tt cz = —02/5 — 3 04/35, 02 = —61 nHz and 04 = —73.5 11Hz. It is thus 
assumed that there is a known negative gradient below the surface down to a 
radius r s = O.95i?0 and the latitudinal dependence of this shear layer is modeled 
by P{9) = (p +Pi cos 4 9)/R Q , so that 



n(r,9) = n c 



Q(9)(r - r cz ) + 5 



[n(r CZl 9)-n c 



SOLA: sdab4.tex; 26 November 2012; 2:20; p. 9 



A. Bonanno 



a 



. 0.6 



0.2 



0.0 



0.8 



0.4 



1.0 



1 .2 



0.6 



0.7 



0.8 




0.2 



04 



Of; 



10 



x 



Figure 4. Left panel: the a-effect profile (solid line) and turbulent diffusivity (dot-dashed 
line) used for an aQ model with a dynamo action produced from subsurface shear instabilities. 
Right panel: the isocontour lines of the differential rotation [f2(r, 6)] in units of the equatorial 
rotation. Note the negative shear near the surface around t/Rq 0.95. 



and Q(6) = (n cq — fl cz + P(9)(Rq — r s ))/(r — r cz ). For actual calculations the 
values po = 437 nHz and p± = —722 nHz have been chosen. In order to confine 
the magnetic field in the subsuface layers it is assumed that % is also maximum 
near the surface, and it sharply decreases by about two orders of magnitude 
just below the supergranulation layer where the a-effect is located. The basic 
features of the (critical) dynamo solution in this case is depicted in Figure [5] for 
a solution with C a = 18, Cn — 1 x 10 3 which implies a period of about four 
years and a turbulent diffusivity r)t = 1.4 x 10 13 cm 2 s~ 1 . The period is clearly 
too small, but this is not surprising as the dynamo wave is basically entirely 
confined in the surface layers as can be seen in Figure The only possibility 
to match the solar period is to further increase the turbulent diffusivity in the 
supergranulation layer, but this would imply that the spatial extension of the 
dynamo wave propagates deeply within the convective zone, thus preventing the 
negative radial shear to produce the correct butterfly diagram. Moreover the 
parity of the solution is clearly symmetric because C a — 15 for quadrupolar 
modes. 

6. Conclusions 

Where is the a-effect located in the Sun? Despite the difficulties present in 
models of flux-dominated dynamo (too low eddy diffusivity, unknown strength 
and location of the return flow), mean-field models with an a-effect located at 




(13) 
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Figure 5. Butterfly diagram for a model with strong radial subsurface shear. The isocontours 
for the toroidal field (upper panel) at r = 0.95R and the radial field on the top boundary are 
shown. Full and dashed lines show the positive and negative levels respectively. The solution 
is obtained with C a = 18 and Cq = 1 x 10 3 which implies a period of about four years and a 
turbulent diffusivity r/t = 1.4 X 10 13 cm 2 s — . 



the bottom of the convective zone are successful in reproducing several aspects 
of the solar activity cycle. In this model the origin of the a-effect is due to 
the turbulent magnetic hclicity and current helicity generated by the quasi- 
interchange instabilities below the tachocline and does not follow from mixing- 
length theory. The most striking result comes from the parity of the solution: 
only with an a-effect located at the bottom of the convective zone the dipolar 
modes are most easily excited, at least at the linear level. This result also apply 
for models with the a-effect generated by the subsurface shear, the solutions are 
mostly symmetric, rather than anti-symmetric. It would be nice to see if a more 
realistic surface boundary c onditio n can solve this problem as proposed by Pipin 
and Kosovichev <|201 lh . 

Although the models discussed in the previous session are kinematic, it is 
difficult to believe that the discussion of the previous session will drastically 
ch ange if a non-linearity via the q-q uenching is included i n the models. Jouve et 
al. ( 20081 ) present a code comparison between kinematic models and non-linear 



models that shows that the critical solutions are basically the same within the 
tolerance of different numerical schemes. On the other hand the advantage of the 
kinematic approach is to provide a complete view of the spectrum of the dynamo 
waves that can be excited for a given set of dynamo numbers. This is clearly an 
important piece of information for the non-linear simulations, for instance for 
investigating degenerate configurations due to the degeneracy in the spectrum. 
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Figure 6. Toroidal and poloidal field configuration for a typical dynamo solution with a 
strong subsurface shear, and C a = 18, Cq = 1 X 10 3 . The left part are the isocontours line of 
the toroidal field, with solid line for negative and dashed line for positive value of the field. 
The right part represents of the streamlines of the poloidal field given by contours of Ar sin 8. 
Solid line are for negative values of A. Note that the dynamo action is confined around the 
region of strongest radial shear, and, for this reason, the period tends to be rather short in 
general (less than ten years in these type of models) . 



An important constraint on the viability of a solar dynamo with a tachoclinc 
a-effect is the fact that successful models must have a positive a-effect: no 
migration is present for a negative a-effect. It would be an important check 
for the theory to show that this is actually the case by means of 3D global 
numerical simulations of the kink and quasi-interchange instability in spherical 
symmetry. Recent investigat ions in this direction in cylindrica l symm etry have 
questioned this possibility ( Gellert. Riidieer and Hollerbachl . l200ll ) although 



further numerical and analytical work is needed before a firm conclusion can 
be reached. 

Acknowledgements I am grateful to Axel Brandenburg and Fabio Del Sordo 
for stimulating discussions. I would also like to thank the organizers of the SDO 
meeting (in particular Sasha Kosovichev) and ESA for financial support. 



SOLA: sdab4.tex; 26 November 2012; 2:20; p. 12 



Solar Dynamo and MHD instabilities 



References 

Bernstein. I.B., Frieman, E.A., Kruskal, M.D., Kulsrud, R.M. 1958, Proc. Roy. Soc. Lond., 
A244, 17 

Braithwaite, J. 2006, Astron. Astrophys. 453, 687 
Bonanno, A. arXiv:1112.4929B 

Bonanno, A.; Elstner, D., Rudiger, G., Belvedere, G. 2002, Astron. Astrophys. 390, 673 

Bonanno, A.; Elstner, D., Belvedere, G. 2006, Astron. Nachr., 327, 680 

Bonanno A., Urpin V. 2008, Astron. Astrophys. 477, 35 

Bonanno A., Urpin V. 2008, Astron. Astrophys. 488, 1 

Bonanno, A.; Urpin, V. 2011, Phys. Rev. E, 84, 056310 

Bonanno, A.; Urpin, V. 2012, Astrophys. J. 747, 137 

Bonanno, A.; Brandenburg, A., Del Sordo, F., Dhrubaditya, Phys. Rev. E 86, 016313, 2012 

and arXiv:1204.0081. 
Brandenburg, A. 2005, Astrophys. J. 625, 539 

Chatterjee, P.; Mitra, D., Rhcinhardt, M., Brandenburg, A. 2011, Astron. Astrophys. 534, 46 
Chatterjee, P.; Nandy, D.; Choudhuri, A. R. 2004, Astron. Astrophys. 427, 1019 
Del Sordo, F., Bonanno, A., Brandenburg, A., Dhrubaditya, M., In: C.H. Mandrini (ed.) Com- 
parative Magnetic Minima: Characterizing quiet times in the Sun and star, IA U Symposium 
286, Cambridge University Press, and arXiv:1111.1742. 
Dikpati, M., Charbonneau, P. 1999, Astrophys. J. 518, 508 
Dikpati, M. and Gilman, P.A. 2001, Astrophys. J. 559, 428 

Dikpati, M., Corbad, T., Thompson, M.J., Gilman, P.A. 2002, Astrophys. J. 575, L41 
Ducz, V., Mathis. S. 2010, Astrophys. J. 517, 58 
Durney, B. R. 2000, Astrophys. J. 528, 486 

Gellert, M., Rudiger G.,Hollerbach, R. 2011, Mon. Not. Roy. Astron. Soc. 414, 2696 
Ferriz-Mas, A. Schmitt, D., Schiissler, M. 1994, Astron. Astrophys. 289, 949 
Gocdbloed, H., Poedts, S., Principles of Magnetohydrodynamics, Cambridge University Press, 
(2004) 

Gruzinov, A.V.; Diamond, P.H. 1994, Phys. Rev. Lett., 72, 1651 

Guerrero, G. and de Gouveia Dal Pino, E. M. 2008, Astron. Astrophys. 485, 267 

Jouve et al. 2008, Astron. Astrophys. 483, 949960 

Kiiker, M., Rudiger, G., Schultz, M. 2001, Astron. Astrophys. 374, 301 

Kitchatinov, L. 2008, Astron. Rep. 52, 247 

Kitchatinov, L., Riidigcr, G. 2008, Astron. Astrophys. 478, 1 

Kosovichev, A.G. 1996, Astrophys. J. 269, 61 

Parker, E. N. 1955, Astrophys. J. 122, 293 

Parker, E. N. 1993, Astrophys. J. 408, 707 

Pipin, V.V., Kosovichev, A.G. 2011, Astrophys. J. 727, L45 

Pipin, V.V., Kosovichev, A.G. 2011, Astrophys. J. 738, 104 

Riidigcr, G., Kitchatinov, L., Elstner, D. 2012, Mon. Not. Roy. Astron. Soc. 425, 2267 and 
arXiv:1107.2548. 

Schou, J., 1991 In: Challenges to Theories of the Structure of Moderate-Mass Stars, Lecture 

Notes Phys. 388, ed D.Gough & Toomre editors (Berlin), Springer, 81 
Silant'ev, N.A. 2000, Astron. Astrophys. 364, 449 
Spiegel, E. A. and Zahn, J. -P. 1992, Astron. Astrophys. 265, 106 
Spruit, H. 1999, Astron. Astrophys. 349, 189 

Steenbeck, M., Krause, F., Radler, K.H. 1966, Z. Naturforsch. 21a, 369 

Tayler R.J. 1973, Mon. Not. Roy. Astron. Soc. 161, 365 

Tayler R.J. 1973, Mon. Not. Roy. Astron. Soc. 163, 77 

Tayler R.J. 1980, Mon. Not. Roy. Astron. Soc. 191, 151 

Zahn, J.-P., Brun. A.S., Mathis, S. 2007, Astron. Astrophys. 474, 145 



SOLA: sdab4.tex; 26 November 2012; 2:20; p. 13 



SOLA: sdab4.tex; 26 November 2012; 2:20; p. 14 



